Spin-Peierls instability in the three-leg Heisenberg ladder 
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Because the three-leg ladder behaves like a renormalized single Heisenberg chain we argue that a 
spin-Peierls instability must occur in this system when it is coupled to three-dimensional phonons. 
Using the bond-mean-field theory, we show that this is indeed the case. The dimerized state below 
the spin-Peierls transition temperature forms into the columnar dimerized phase not the staggered 
one. This contrasts with the argument based on antiferromagnetism. A physical argument based 
rather on spin bonding into singlets explains why the columnar configuration is favored. No quantum 
criticality (gaplessness) can occur in the columnar arrangement of the dimerized chains. 
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Interest in the Heisenberg ladders stems from one of 
their most intriguing aspects: namely, the critical prop- 
erties are dependent on the number of legs, and also from 
the various (or at least potential) exotic quantum states 
that can be realized in them (for a review see Refs. [HQ. 
The spin excitations in the n-leg ladders with an even 
number n of legs are gapped, while those with an odd 
n have gapless excitation spectra The even- n-leg 
ladders provide an example of spin liquids. The latter, 
defined (loosely) as gapped phases without long-ranged 
order (LRO), are, in particular, believed to be relevant 
to the physics of high-T c superconductivity Q . 

In this Letter we study the effect of phonons on the 
three-leg ladder. Having an odd number of legs, this 
sytem is gapless in the absence of coupling to phonons. 
This was confirmed, e.g., by the Monte Carlo simulations 
[1, [j| . In the limit of strong rung coupling the three- leg 
ladder behaves as an effective (renormalized) Heisenberg 
chain [l[ . A recent thorough theoretical study on this sys- 
tem corroborated this point [6J . Experiments on the real 
three-leg ladder compound Sr 2 Cu305 find its spin suscep- 
tibility similar to that of a (gapless) Heisenberg chain [7( . 
It is therefore natural to predict that a spin-Peierls (SP) 
instability can take place in the three-leg ladder when 
it is coupled to phonons. When a spin- 1/2 Heisenberg 
chain is coupled to phonons, lattice distortions and a spin 
gap occur simultaneously below a transition temperature 
Tsp known as the SP temperature Q. The gap is due 
to the distortion-induced dimerization. The SP systems 
have been extensively studied both experimentally and 
theoretically @, 0, [13, EH - The SP instability has been 
observed in TTF - CuBDT^ TTF - AuBDT [H, [H, in 
other organic compounds [101 ] , and in the inorganic com- 
pound CuGe0 3 [14 1 . 

We use the Jordan- Wigner (JW) transformation [l5| 
and bond mean-field theory (BMFT) [HElEl to show 
that indeed a SP transition takes place in the three-leg 
ladder when it is coupled to three-dimensional phonons. 
To the best of our knowledge, there are no other theoreti- 
cal studies of this problem, and real three-leg ladders with 
an SP instability do not exisit yet. There are, however, 



studies of the intrinsically dimerized ladders. According 
to Refs. 



18 



19j . for the antiferromagnetic three- leg lad- 
der with dimerization preset in the staggered pattern (cf. 
Fig. Hk), there exist a critical line in the dimcrization- 
rung-coupling plane (<5, J_l), where the system is gapless. 
Dimerized ladders provide a counterintuitive example of 
"restored quantum criticality" , when a system (ladder) 
built from gapped blocks (dimerized chains) can be gap- 
less, contrary to naive expectations. See Refs. [HI, [l§] 
and Refs. therein for more detail. Our analysis shows 
that this interesting phenomenon does not occur in the 
coupled spin-phonon ladder: when the ladder is allowed 
to choose the dimerization pattern from the minimum 
energy condition, it orders into the columnar phase (cf. 
Fig. [lb)- In the latter case, the ladder of dimerized chains 
is always gapped ■ So, the SP transition in a three- leg 
ladder qualitatively resembles that in a single chain. 

The spin Hamiltonian for the three-leg ladder with an- 
tiferromagnetic couplings and coupled to phonons is 
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where i is the site label along the chains (i.e., rungs), 
j = 1, 2, 3 labels the legs (chains), and N is the number 
of sites in a single chain; the total number of spins is N t — 
3N. J± is the coupling along the rungs and Ji,i+i(j) ^ s 
the longitudinal position-dependent coupling because the 
chains of the ladder are linearly coupled to the phonon 
field u^i+i^') 2]. Restricting our analysis to the static 
alternating lattice deformations u along the chains, we 
take 

Ji.t+iti) = Jo{l + l(ui,i+i(j))), (« M+ i(j)) oc {-l) l u 

In this (adiabatic) approximation the phonon Hamilto- 
nian is given by its static deformation part 

1 N t J S 2 



H ph = -N t Ku^- — 



(3) 



where the dimensionless dimerization parameter S = ju 
and the spin-phonon coupling A = Jq7 2 / K. The alter- 
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FIG. 1: The bold/thin solid lines represent the 

stronger /weaker chain coupling Jo(l ± <5), respectively. The 
dashed lines correspond to the rung coupling J±. The dimer- 
ization pattern is staggered in (a) and columnar in (b). 



ary conditions along the rungs, and using the Nambu 
formalism, the BMFT yields the single-particle effective 
Hamiltonian 



(6) 



where the Hamiltonian density H^ co ^ is a 6 x 6 matrix 
given by 



nated frozen intrachain displacements @ result in dimer- 
ization of each of the three chains. The dimerization of 
the whole ladder can be in the staggered or columnar 
patterns, as shown in Fig. [1] 

The total Hamiltonian assumes then the form H = 
Hsl + H p h, and the effective intrachain spin coupling 
dependent on the dimerization pattern is 



Ji, 
Ji, 



Jo[l + (- 
Ml + (- 



-iy +J 5] staggered 
-iy5] columnar 



(4) 
(5) 



The BMFT we apply to the JW fermions consists of 
two key approximations. First, in dealing with the phases 
resulting from the JW transformation, the phase differ- 
ences due to hopping of the JW fermions around any 
given elementary plaquette is set to be equal to n. Sec- 
ond, the quartic JW fermionic terms c\ jCijc\ +1 jCi+x ,j 
resulting from the Ising interactions are decoupled using 
the bond parameters [y, 11, ID, 1(1, 17|. The latter are 
defined as Q+ = (c2i,j(^ i+ltj ), Q- = (c2i+i,j4i+2,j) for 

j = 1,3; and Q' + = (c 2 i )2 4i+i,2) > Q- = (^21+1,2^+2,2) 
for j = 2. For chains 1 and 3, the same bond param- 
eters are used because the ladder is symmetric under 
exchanging chain labels 1 and 3 p\. In the direction 
along the rungs, only one bond parameter is sufficient, 
P = ( c 2i,jcli,j+i)- For the reasons to be explained be- 
low, we will concentrate on the columnar dimerization 
pattern ([5]) . Fourier transforming along the chains direc- 
tion, keeping the chains labels because of the open bound- 
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Here c" fc is the Fourier 
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transform of c°j along the chain j; i.e., with respect to 
the index i (a = A, or B), and because of the antifer- 
romagnetic correlations the lattice is subdivided into 
two sublattices A and B. The parameters entering the 
effective Hamiltonian are: A* = (j| + e ifc - jjLe" lfc )/2, 
and C = Jii/2, where 

Jf ± = J (l±5)(l + 2Ql) , J ±1 = J±(1 + 2P), 



C t = NJ + \Q + \ 2 + NJ_\Q_\ 2 + ^-J_\Q 



I |2 



^J+\Q> + ? + 2J+N\P\ 2 + i -^f. (8) 



Here, j\ ± = J%± for Q± = Q or J' 1± for Q± = Q'. Diag- 

onalizing H^ co ' yields six energy eigenvalues ±Ej CO \k), 
j = 1, 2, 3, where 



(9) 



and 



J 
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n = 2,3. 



(10) 



In Eqs. (|9TT0f 

z> = (j[ + f + {,][_) 2 -2j[ + j[_cos{2k), 
t = (Ji+J{_ + Ji-J[+) cos(2fc) - (Ji+J' u 



Ji-J[-). 



In the absence of dimerization (8 = 0), the energy eigen- 
values (|9I10[) coincide with those found in Ref. [6j. 
The partition function of the single-particle Hamilto- 



nian (JGJ) can be calculated and leads to the following free 
energy per spin 

3 



F (co) = Q_ 
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where (3 = 1/fcsT. Finally, the mean-field equations are 
derived from minimization of the free energy (llip with 
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FIG. 2: The free energies of the three-leg ladder without 
phonons (A = 0, no dimerization) , and in the presence of 
phonons (A = 0.1) for the staggered and columnar configura- 
tions are plotted as functions of ftemperature. 



respect to the six mean- field parameters Q± , Q ± , P and 
S. These self-consistent (integral) equations are solved 
numerically. 

The mean-field equations predict no magnetic long 
range order (LRO), even at zero temperature. At the 
same time, they predict the simultaneous appearance at 
some critical temperature of the structural LRO (lattice 
dimerization 5) and the spin gap, generated by 5 =/= 0. 
So, this is a SP transition. 

First, we compare the free energies of the two dimeriza- 
tion patterns. For the staggered configuration with the 
effective coupling (|4|) , the single-particle effective Hamil- 
tonian, its spectrum, and the mean-field equations are 
derived in the same manner as described above for the 
columnar configuration. Fig. fj] shows the free energies as 
functions of temperature for both dimerized configura- 
tions and without dimerization (A — 0). The free energy 
of the dimerized columnar configuration is lower than 
those of the staggered or non-dimerized configurations 
below a nonzero temperature, which we identify as the 
SP temperature T$p- So, the columnar order constitutes 
the thermodynamically stable state at low temperature. 
Above T$p the dimerization disappears, and all three free 
energies become equal. We also plot 
as a function of a 
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ture in Fig. [3] for A = 0.5. 

passes through a maximum, then decreases rapidly in the 
strong coupling regime a»l. This means that the stag- 
gered and columnar states become practically degenerate 
when a»l. 

To understand why the columnar configuration is the 
stable one, assume a is large. In this limit, the three-leg 
ladder behaves pretty well like a single Heiscnberg chain 
with an effective coupling, and the spin degrees of free- 
dom on two of the chains freeze into spin singlets on the 
rungs @. It is then favorable for the coupling along the 
chains to couple the singlets into a plaquette (which gives 
rise to a linear combination of transverse and longitudi- 
nal spin singlets for the spins on four of the six sites) 




FIG. 3: The free energy difference between the staggered and 
columnar configurations as a function of a = J±/Jo at T = 
0.009 Jo, A = 0.5. 



0.4 
0.3 
0.2 
0.1 



0, 



— V 1 1 1 1 
■ \ (a) 


1 I 1 




— e + " 


- ' 


— e 


/ \ 








i . i 


i 







0.5 1 1.5 

k,.T/J n 




FIG. 4: The temperature dependence of the mean-field pa- 
rameters Q±, Q'±, P and 5 for a = 1, A = 1. 



like in Fig. QJb). The staggered configuration would be 
favorable if the spins on the rungs were ordered antifer- 
romagnetically. For intermediate and weak values of the 
rung coupling a, the spins are not frosen but are nonethe- 
less locked into random singlets on two of the three chains 
@, and the argument based on spins singlets still holds, 
while that based on antiferromagnetically ordered spins 
on the rungs continues to be wrong. 

The mean-field parameters are plotted as functions of 
temperature in Fig. [4] for A = 1 and a = 1. At T$p, a 
spin gap opens and a structural phase transition occurs 
from the disordered phase into the dimerized phase. This 
is clearly indicated by the temperature dependence of S 
in the inset of figure [4] (b) , which is very similar to that 
in an ordinary second-order phase transition. The other 
(bond) parameters are not critical. Below Tsp, Q\ ^ Q! 
due to non-zero dimerization, while = QL above T$p. 

The SP instability is accompanied by the opening of 
an energy gap in the spin excitation spectra (spin gap). 
The spin excitation energies consist of three bands Ej (fc), 
Eqs. (fT0| . with E 2 and E 3 gapped even in the ab- 
sence of dimerization with a gap of the order of J± . The 
gap in E\{k) is induced by the dimerization i5. We plot 
E\(k) for the columnar dimerized (A = 0.7) and non- 
dimerized (A = 0) ladder in Fig. [5l The spectrum of the 
non-dimerized three-leg ladder, calculated in Ref. Q is 
gapless, i.e., Ei(Q/tt) — 0, as it should be for an odd-leg 
ladder (Note the symmetry of Ej(k) with respect to 
the point k — The spin gap E g = Ei(0) defined from 
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FIG. 5: The spectra for the uniform (A = 0) and columnar 
dimerized (A = 0.7) ladders are shown. Here, a = 1. 



state can occur for such arrangement of the dimerized 
chains. Finally, we conjecture that the results derived 
here for the three-leg ladder should stay true for any n- 
leg ladder with n = 5, 7, ... 
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Eqs. 



E, 



and © is 



Q+ + Q-) + (Q+-Q-)]. (12) 



So, to leading order our theory, as any mean-field one 
[2 11 ] gives the non-interacting fermion (equivalcntly, the 
dimerized single XU-chain) result E g oc 5, albeit renor- 
malized by the bond parameters Q±. This is confirmed 
by direct numerical calculations of E g (8). Note that 
for a dimerized chain E g oc 5 2 ^ 3 [E|, c orrected by the 
marginal logarithmic prefactor |lOl 23. [23| . The tem- 
perature dependence of the spin gap is shown in Fig. [5] 
Both the numerical calculations and qualitative analyses 
of the mean-field equations show that <5, E g , and T$p 
increase monotonously with the phonon coupling con- 
stant A. However the dependence of those parameters 
on a = J±/Jo is trickier. The parameters 5, E g and Tsp 
initially decrease with growing a, but then saturate as 
a > 5 as seen in Fig. [S] The mechanism of such satura- 
tion was already discussed [l|, [|| , and can be understood 
in the regime q > 1 as due to the fact that the ladder 
behaves as a "renormalized" spin-^ chain. 

In this work, the spin-Peierls instability in the three- 
leg ladder coupled to phonons is studied. We map the 
spin operators of the three-leg Heisenberg ladder onto 
Jordan- Wigner fermions. The resulting interaction terms 
are decoupled within the bond-mean-field theory. This 
theory yields a spin-Peierls transition into the columnar 
dimerized phase. Qualitatively, the three-leg dimerized 
ladder behaves like a single spin-Peierls chain. No gapless 



0.06 







Q. 




GA 


0.04 > 


> 




O 




cr 


0.02- 


LU 




z 




LU 





Dn °n 



— 1 ' 1 



oc= 



° 1 

a □ a 5 
°o ft °30 



a 
2 



0.0 



0.06 



k B T/J 



FIG. 6: The temperature dependence of the spin gap is shown 
for several values of a, and for A = 0.1. 
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